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Abstract. The recent results on various Beyond the Standard Model (BSM) Higgs boson searches 
performed by the D0 experiment at the Tevatron are presented here. In particular, the Higgs bosons 
in supersymmetric models and fermiophobic scenario have been investigated. No significant excess 
over the Standard Model (SM) expectations have been observed and accordingly limits have been 
established on the corresponding model parameters. 
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INTRODUCTION 

Besides the Standard Model Higgs Boson searches, the D0 physics program is enriched 
with various well motivated BSM Higgs boson searches. The searches are primarily fo- 
cused on the Higgs bosons from three different models: the Fermiophobic Higgs model, 
the Minimal Supersymmetric Standard Model or MSSM and the next to Minimal Super- 
symmetric Standard Model or NMSSM. The analyses summarized here are performed 
with 1-4.2 fb 1 of pp collision data recorded with the D0 detector during Run II of the 
Tevatron Collider at Fermilab. 

FERMIOPHOBIC HIGGS SEARCH 

The Fermiophobic model assumes zero couplings of the Higgs boson to fermions while 
the Higgs couplings to the gauge bosons remain the same as in the SM. In such models, 
Higgs boson production via gluon fusion is absent and the Higgs bosons are mainly 
produced in association with a W ± or a Z° boson (VH) and via vector boson fusion 
(VBF). D0 has performed the searches for such Fermiophobic Higgs bosons (hf) in two 
different channels: hf — > 77 and W ± hf(^ W + W~) with 4.2 fb -1 and 3.6 fb -1 datasets 
respectively. 

In the inclusive hf — > 77 search, both both VH and VBF production modes are 
considered. The background events for this channel are contributed by Z°/t* — > e + e~, 
74-jet and QCD dijet production processes. Here the invariant mass (MJ distribution 
from two reconstructed photons (with px > 20 GeV) has been utilized as the final 
discriminating variable. The observed limits on the production cross section have been 
translated into the benchmark model expectations and a Fermiophobic Higgs boson mass 
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FIGURE 1. (a) The limits on Br(hf — > 77) as a function of Higgs boson mass (Mh f ); Mh f < 102.5 is 
excluded at 95% CL. (b) The 95% CL limits on a(pp ->■ W ± h f ) x Br(h f -> W+W") as a function of 
Higgs boson mass. The previous D0 results are also shown for comparison. 



below 102.5 GeV is excluded at 95% CL (see Fig. (Ha)). 

In the Fermiophobic Higgs models, the BR(hf — > W + W~) becomes nearly 100% 
for a Higgs boson mass above 100 GeV [Q]]. The signatures of W ± hf(^ W + W~) 
events have been looked for in the events containing two like- signed leptons, 11 
(electrons or muons). The physics processes like W ± Z° — >7 / and Z°Z° —titbit* 
along with the misreconstructed Z°/y* — > e + e~//i+/i~ events are the most dominant 
background processes for this search channel. The observed and expected limits on 
c(pp — ► W ± hf) x Br(hf — ► W + W~) are compared to the SM and Fermiophobic model 
predictions in Fig.QJb). 



SEARCH FOR NEUTRAL MSSM HIGGS BOSONS 

The MSSM^lOl] requires two doublet Higgs fields to generate masses to both "up"- 
type and "down"-type fermions. After electroweak symmetry breaking, such a two- 
Higgs-doublet model predicts the existence of five physical Higgs bosons: two CP- 
even Higgs bosons, h and H, one CP-odd Higgs boson, A and a pair of charged Higgs 
bosons, H ± . At the tree level, the MSSM Higgs phenomenology can fully be described 
by two free parameters: niA, the mass of the CP-odd Higgs boson and tan/3, the ratio 
of vacuum expectation values of the Higgs fields to the "up"-type and "down"-type 
fermions respectively. 

The couplings of neutral Higgs bosons to bottom quark (b) and tau lepton (t) i.e., 
"down"-type fermions scale as tan/3 with respect to their SM values. Therefore for high 
tan/3, the production cross sections involving b quarks are enhanced of a factor tan 2 /3. 
Moreover, the CP-odd Higgs boson A becomes degenerate with either one of the other 
neutral Higgs bosons, h/H which leads to a total 2 • tan 2 /3 enhancement in production 
cross section relative to the SM prediction. At high tan/3, the neutral Higgs bosons 
(h/H/A = 4>) decay predominantly into bb (Br^90%) or t + t~ (Br^l0%) pairs. 

Three different analyses viz., d> — >• T + T~, 4>b(b) — > T + T~b(b) and d>b(b) — > bbb(b) 
are presented here consisting of 2.2, 1.2 and 2.6 fb -1 datasets respectively. The inclusive 
4> — > t + t~ search considers the events where one of the tau leptons decays into a muon 
(Tjn) while the other tau lepton decays hadronically (ihad)- The hadronically decaying tau 
leptons are classified into three different types depending on their detector signatures and 
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FIGURE 2. (a) The visible mass distribution (see text) for inclusive pp — > <!>(—» t + t~) searches, 
(b) Limits on cr(pp — * 4>b(b)) x BR(4> — > t + t~) vs M<j>. (c) The tree level limits on tan/3 — hia from 
<t>b(b) — > bbb(b) search channel. 



are selected through the usage of neural network discriminants (NN T ). The background 
events for this channel are dominated by Z°/y* — > t + t~, Z°/y* — > and multijet 

QCD production processes. The W ± +jet events are rejected by requiring Mt(Pj,Et) < 
40 GeV. The said analysis with 1.2 fb 1 Run lib dataset is combined with earlier 
analysis 0] with 1 fb -1 of Run Ha dataset which considers additional final states viz., 
T e Thad an d TeT^. The distribution for the combined visible mass i.e., a/ (PtI + Pt2 + Pt) 2 > 
where P T i 2 is the four vector of the visible tau decay products and Pp = (Bp, Ex, My, 0), 
is shown in Fig. Ufa). In the 4>b(b) — ► T + T~b(b) search, the events are characterized 
by an isolated muon, a hadronic tau jet candidate (same NN T selection as in 4> — > 
T + T _ ) and a b-tagged jet. Multijet QCD, tt and Z°/y* — > T + T~b/c are the dominating 
background processes and they are reduced further by applying a 2D neural network 
discriminant. Fig. |2£b) shows the limits on a(pp — > <3>b(b)) x BR(4> — > t + t _ ) at 95% 
CL. In the 4>b(b) — > bbb(b) search, the events are required to have 3-5 jets (px > 20 
GeV and \rj\ < 2.5) three of which are b-tagged. The overwhelmingly large QCD 
background contributions are determined through usage of both data and Monte Carlo 
(MC) distributions. The invariant mass distribution from the jet pairs are considered 
here after applying a likelihood discriminant for further separation between signal and 
background processes. The cross section limits are translated into the tan/3 — niA space 
assuming a simple 2 ■ tan 2 /3 enhancement over the SM and are shown in Fig [2c). 




FIGURE 3. Combined limits on MSSM parameter space for different benchmark scenario: 
m™ x ,jU=-200 GeV (left) and no-mixing, n=+200 GeV(right). 

Finally, all three above mentioned searches are combined taking into account system- 
atic uncertainties, including the correlations across channels and analyses. Fig. [3] shows 
the combined results interpreted in different benchmark scenarios oD. 



NMSSM HIGGS SEARCH 



In the next-to-MSSM (NMSSM) [6], at lower masses the Higgs boson predominantly 
decays into a pair of lighter neutral pseudoscalars, a. The result is a suppression of h — ► 
bb branching ratio. Depending on its mass, the neutral pseudoscalar can decay into a pair 
of muons (2M M < M a < 3M K ) or taus (2M T < M a < 2Mb). Thus the final state signatures 
would consist of two muon pairs or two tau pairs. However, because of the difficulties 
with the 4t final state, the search is performed in the h — > a(— > T+T~)a(— > /i + /i~) final 
state. Due to the extreme collinearity of the pair from the pseudoscalar decay 

and the finite angular resolution of the D0 muon system, each muon pair in 4 muon 
channel is required to contain one muon trajectory associated with a companion track. 
A requirement of AR > 1 between the muon pairs is applied to ensure that these muon 
pairs are well separated from each other. For the 2t + 2fi channel however both the 
muons in a muon pair are required to be reconstructed and should be close to each 
other within AR < 1. It is to be noted here that 2r + 2/1 analysis does not require 
explicit tau reconstruction (compromised by the collinearity of the taus). Instead, the 
Ex requirements are used to predominantly select the events with tau candidates. 

For M a < 2M T , the 4 muon channel sets an upper limit of 10 fb on c(pp — > hX)Br(h — > 
aa)Br(a -> /i + /i~) 2 . Assuming Br(h -> aa) ^100% and the SM cr(pp -> hX) « 1.9 pb 
(for Mh = 100 GeV), the corresponding upper limit on Br(a — > /i + /i~) « 7%. However, 
in the NMSSM the BR(a — ► /i + /i~) is expected to be larger than 10% for M a < 2m c 
and depending on the branching ratio of a — > cc possibly even up to 2m T . Thus the 
recently reported[7] D0 result severely constrains the region 2M jti < M a < 2M T . For 
M a > 2m T , the limits set by the 2t + 2jU channel are a factor of 1-4 larger than the 
expected production cross section. 

CONCLUSIONS 

D0 has made rigorous BSM Higgs bosons searches in most of the promising channels. 
Neither of the searches presented here observe a statistically significant excess over the 
SM background expectations. Accordingly at 95% CL, the most restrictive limits have 
been established on relevant parameters for different theories. Improved analysis tech- 
niques along with the rapidly increasing datasets are expected to enhance the sensitivity 
of future BSM Higgs searches. For the most recent results see Ref. Jsj]. 
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